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Abstract 

Particle production in high-energy collisions is often addressed within the frame- 
work of the thermal (statistical) model. We present a method to calculate the canon- 
ical partition function for the hadron resonance gas with exact conservation of the 
baryon number, strangeness, electric charge, charmness and bottomness. We derive 
an analytical expression for the partition function which is represented as series of 
Bessel functions. Our results can be used directly to analyze particle production 
yields in elementary and in heavy ion collisions. We also quantify the importance 
of quantum statistics in the calculations of the light particle multiplicities in the 
canonical thermal model of the hadron resonance gas. 



1 Introduction 



The statistical model of the hadron resonance gas was shown to be very successful in ad- 
dressing particle production in high-energy collisions [1.2,3.4.5.6.7.8.9,10. 11, 12J. In central 
nucleus-nucleus collisions the formulation of conservation of quantum numbers (sometimes 
genericaliy referred to as charges) in the grand canonical (GC) ensemble is well suited 
[2l[3l4l5ll6] . However, for small systems like e + e~ or pp as well as for peripheral nucleus- 
nucleus collisions only the canonical (C) ensemble gives the correct description of hadron 
yields pTf7]l8l9irT0|lllll2j . In the GC formulation the conservation of charges is implemented 
on the average and is controlled by the appropriate chemical potentials whereas in the 
C-ensemble charges are conserved exactly [1J. 

The exact treatment of quantum numbers in statistical mechanics is well established 
|13fl4|ll5|16|ll7fl8|19j . It is in general obtained by projection of the partition function on 



the desired values of the conserved charges by using group theoretical methods [T5lll7|ll8f20ll21 
For studies of hadron production one needs in general to account for the exact conser- 
vation of five quantum numbers: the baryon number N, strangeness S, electric charge 
Q, charmness C, and bottomness B. The conservation is a direct consequence of the 
G — C/jv(1) <8> <8> Ub(1) invariance of strong interactions. 

The implementation of an exact conservation of quantum numbers in the canonical parti- 
tion function for the hadron resonance gas requires integration over the symmetry group 
[Tlfl5] . This in general leads to numerical problems due to the oscillatory behavior of the 
integrand. To avoid numerical problems, methods were developed to express the partition 
function as series of Bessel functions which is well suited for numerical implementations 
[2"2"|23f24ll2"5] . An explicit analytical expression for the canonical partition function using 
the above methods was obtained up to now only for three conserved quantum numbers 
i.e. for N, S and Q and employing Boltzmann statistics [24]. 

In the application of the statistical model to hadron production, e.g. in e + e~ annihilation, 
one must employ the formulation of the canonical partition function which accounts for 
the conservation of all five quantum numbers including charmness and botomness since 
hadrons carrying heavy quarks contribute sizably to the overall hadron yields [71l^lllUf26] . 
In addition, for e + e~ collisions, the data are measured with high precision; see for 
a recent summary. Consequently, in the statistical model analysis the approximation of 
particle momentum distributions using Boltzmann statistics is not sufficient [91127] . 

The main scope of this paper is the extension of the previous results for the canonical 
partition function for a hadron gas [23] to include the exact conservation of all quan- 
tum numbers carried by the light and heavy quarks [28]. We will formulate the canonical 
partition function as series of the modified Bessel functions. Our expression is numeri- 
cally stable and can be used to quantify thermodynamics independently of the values of 
the thermal parameters and the initial quantum numbers. Our result accounts for the 
quantum statistics for bosons. For all fermions the masses are significantly larger than 
the temperature. Consequently, the implementation of Boltzmann statistics for fermions 
is a very good approximation. An extension of our results to quantum statistics for all 
particles including fermions would however be quite straightforward. 

As model implementation of the canonical partition function we illustrate its applica- 
tion to the calculations of hadron yields. We quantify the quantum statistics effects for 
the multiplicity of pions and kaons. We also discuss deviations from the exact quantum 
statistics results when the Boltzmann approximation is used for different values of the 
thermal parameters. The hadron resonance gas we are working with is containing all 
known hadrons, including the multi-strange hyperons up to S = ±3 and all charmed and 
bottom hadrons, as listed in the latest compilation by the Particle Data Group ]26j. 

The paper is organized as follows: in the next Section we derive the analytical expression 
for the canonical partition function that accounts for the conservation of five quantum 
numbers. In Section 3 we present a particular numerical implementation to calculate the 
pion and kaon multiplicities and discuss the role of quantum statistics. We summarize our 
results in Section 4. 



2 The canonical partition function 



The appropriate tool to deal in a statistical mechanics framework with a system of quan- 
tum numbers X = (N, S,Q,C, B) related with the G = C7jv(l) ® .. <g> Ub(1) symmetry 
group is the canonical partition function [T7j (h — c — 1): 



■ixj<f>\±l 



where the vector X=(N, S, Q, C, B) characterizes the initial quantum numbers of a system 
related with each of the U(l) symmetry groups and the ± sign refers to fermions and 
bosons. The exact quantum number conservation is implemented by the integration over 
the group G with <fi = (<p N , <pg, 4>Qi 4>c, 4>b) being an element of G. The vector Xj = 
(Nj, Sj, Qj, Cj, Bj) describes the quantum numbers of a particle j with mass m,j and the 
spin-isospin degeneracy factor gj. The sum in the exponential is taken over all particles 
and resonances which carry quantum numbers related with an internal symmetry G. 

The partition function (TfJ depends only on two parameters, the temperature T and the 
volume V of the system. However, the integral representation of this partition function 
is very inconvenient for numerical implementations. This is particularly the case for non- 
vanishing initial quantum numbers due to the oscillatory nature of the integrand. An 
additional complication appears due to quantum statistics effects. Clearly, one could sim- 
plify the problem and employ the Boltzmann approximation by keeping only the first 
term of the series expansion of the ln(l ± x) ±x in the exponential. However, such an ap- 
proximation is only valid if the particle mass is significantly larger than the temperature. 
As the values of temperatures extracted from the fits of hadron abundances obtained in 
heavy ion and elementary collisions [Tf2TITf9l[T0] are close to the pion mass, the Boltzmann 
approximation will give rise to deviations from the correct quantum statistics values in 
particular for light particles like pions or kaons. To calculate the partition function accu- 
rately (using quantum statistics) one needs in general to include the whole series of the 
expansion of the logarithm 

oo k 

ln(l±x) ±1 = 5:(±l)* +1 ^- (2) 

k=i K 

in the partition function for all particle species. However, as the fermions are heavy, the 
error caused by the Boltzmann approximation in this case is very small. Therefore we 
will use the whole series of the logarithm only for light bosons (mesons without charm 
or bottom quarks) and apply the usual Boltzmann approximation to fermions and heavy 
bosons (bosons containing charm or bottom quarks). With this approximation equation 
(CQ) becomes: 

Z n , S ,q,c,b{X) = T^yJ ^ 5 0e^exp|E4 e "" J ^ + EE^^4 (3) 



b k=2 



with the particle partition function 



/ 7 y2_i_ tt1 2 



where j runs over all particles, while b runs only over the light bosons. 

For charm and bottom hadrons, because of their large masses the term y\- z) is much 

. . 7 ° ^'Jc/b Jc/b 

lesqjj than 1. Consequently, we can use the following approximation [7]: 

(.Jc/b ) jc/b 

Inserting now Eq. ([5]) into Eq. ([3]) we obtain: 

Znaq.cAX) « t^tj / rf 5 e"*e«*) 

jc b & l Z7F J J 

with 

oo 

f{$) = £ 4 e ~"^ + EE afc-«** (7) 

j b k=2 

The index j in Eq. ([?]) runs over all hadrons except those which carry heavy flavors 
whereas in Eq. ([6]) j c runs over all charm and j& over all bottom hadrons. 

From Eq.([6]) it is transparent that the integrals over <pc and (f>B related with the charm 
and bottom quantum numbers contribute to the partition function as Kronecker delta 
functions. Consequently, 



+ E E 44(2^ / ft^^ettoSc^SB^+B. 



jc ib & 
c Jb=° 



1 The particle partition function zj ^ is O(10 4 ) for charm particles and 0(10 13 ) for bottom 
particles 



where X and Xj are now three-dimensional vectors composed of the baryon number, the 
strangeness and the electric charge, while the charmness C and bottomness B appear 
only through the Kronecker functions. 

With the approximation dSJ) for charm and bottom contributions to the generating func- 
tional we have reduced the five dimensional integrations to three dimensional ones in the 
canonical partition function. The integrals to compute are of the following generic form: 

lN,S,Q = 7^ r d ^ ( 9 ) 

{iny Jo 

and correspond to the canonical partition function with the conservation of three quantum 
numbers that accounts for the quantum statistics of bosons. 

Table 1 

The combinations of quantum numbers (iV=baryon number, S=strangeness, (5=electric charge) 
for hadrons and the corresponding notation for the sum of particle partition functions for each 
hadron class, Z^adr- The correspondence of Zhadr to the index rij in Eq. [13] (see text) is also 
given. 

Quantum numbers Z^adr Index rij 
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z 
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Z K o 
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s=o 
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Z TT± 
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Z p 
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Z K ± 
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n 5 
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Z E - 


n 7 
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S=-2 


Q=0 


Z S o 


n 8 


N=l 


S=-2 


Q=-l 


z E ^ 


n 9 


N=l 


S=-3 


Q=-l 




"-10 



The integral representation of In,s,q is n °t convenient for numerical analysis as the in- 
tegrand is a strongly oscillatory function, particularly for large initial quantum numbers 
N, S or Q. Following the methods described in Refs. [221123] and [21] we express Tn,s,q 
as a series of Bessel-functions. First, we observe that, in the argument of the exponential 
function in Eq. ([HD, the particles appear pairwise with their anti-particles. Second, the 
contributions of all particles in the sum can be grouped into 14 categories defined by their 
quantum numbers, see Tabled] for definitions. For instance Z^a describes the sum over all 
partition functions z] of particles j with N=0, S=l and Q=0. Consequently, we rewrite 



the integral dHJ) as follows: 



1 /" 27r 

1 n ^q = exp(Zo)— / d(j> N e lN ^ exp[Z r 
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iS(t>s exp[Z K0 (e i(t>s + e-^ 3 )] 
— I . i( ,Qe lQ ^exp[Z w± (e^ + e- li 

Z7T JO 

exp[Z p (e i( ^ + ^ ) + e ^ ( ^+^ ) )] 
exp[Z AT (e i(0JV ~^ ) + e ^(^-^Q))] 
exp[Z A++ (e i( ^ +2 ^ ) + e - i (^+ 2 ^))] 
exp[Z^±(e i(</,s+ ^ ) + e -W*+*a))] 

cv„[7. ( J(<t>N-<t>s) i -i(<t>N-<t>s)\] 



exp[Z K ±{e' . M 

exp[Z A (e i(</,JV - 0s) + e _i( ^-^)] 
exp[Z E+ (e i(fe - </,s+ ^ ) + e-Ww-^+to))] 
exp[Z s -(e i(< ^ JV - </,s -^ ) + e -*(^-^-* ))] 
exp [Z^o 
exp[Z STv - 

r oo 
.k=2 



exp[Z s -(e i( ^-^-^ 
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exp 
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(10) 



Applying the relation: 



exp 









K)] 


.2 





E * m/ « 



(11) 



and the integral representation of the Bessel-function of order h: 

1 r 2 * 

Ih(x) — tt" / rf0exp(xcos0) exp(— i/i0), 

27T JO 

the group integrals in Eq. (JTU1) can be done explicitly, yielding: 



(12) 



10 



1n,s,q = exp(Zo) ]J 

3=1 



E r^pz, 



hadr . 



Tii=—oo 
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fc=2 



E i™. k {2z*±) 



m^=— oo 



n 

h=2 



E -^(^x 



rrih=— oo 



n 
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mi=—oo 
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(13) 



with 



v \ — Q + n i ~ n 2 + 2n 3 + n 4 + n 6 — n 7 — n 9 — n 1Q + 2J ^w-fc + 
1^2 = 5 + ri4 — 715 — n6 — n-j — 2n§ — 2ng — 3nio + ^ hrrih + Zm; 



f 3 = N + rii + n 2 + n 3 + n 5 + n 6 + n 7 + n$ + n 9 + nio, 

and with the specific pertaining to each index rij according to Table [TJ The above 

equation is a generalization of that for a canonical partition function with exact conser- 
vation of baryon number, electric charge and strangeness [24J to the case where quantum 
statistics for bosons is explicitly included. 

To obtain the partition function conserving five quantum numbers we use Eq. ((9]) to 
transform Eq. (jSj) into: 

Zn,s,q,c,b(X) ~%n,s,q 5c,o5b,o 



This expression, together with (|13p is our final result for the partition function that 
accounts for exact conservation of baryon number, electric charge, strangeness, charmness 
and bottomness. The partition function (JHJ), contrary to its integral representation (PQ), 
is free from oscillations and is numerically stable independent of the values of the thermal 
parameters or the values of initial quantum numbers. 

The partition function (|14p can be used to describe thermodynamical properties of the 
hadron resonance gas under constraints of the exact conservation of all relevant quantum 
numbers. In particular, from Eq. (fT4"j) we obtain the multiplicity (rij) for hadron species 
j by introducing a fugacity parameter Aj which multiplies the particle partition function 
Zj and by differentiating: 



(14) 



+ E E Z i Z j b X N+N JC +N Jb ,S+S JC +S n ,Q+Q ]c +Q Jb Sc,C ]c SB,B ]c +B : 





dlnZ 



■>N,S,Q,C,B 



(15) 



For bosons, e.g. for n^, 



one obtains from Eqs. (fT4"l) and ffT5l) . 



(«W±) 



dlnZ 

<9A„-± 



7T- 1 - 

Zn,s,q,c,b{X - kx^± 

t =; 

Zn,s,q,c,b{X) 
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(16) 



k=l 



whereas for fermions there is only one term contributing (k—1) because fermions are well 
approximated by Boltzmann statistics. 



3 Numerical results 



We have applied the above obtained analytical expression for the canonical partition 
function (|14p to quantify particle production in e + e~ annihilations at LEP energies [9f29] . 
We will not repeat the discussion given there but focus, in the following on an illustration 
of the importance of quantum statistics in the calculation of multiplicities of light bosons. 




Figure 1. The deviation of the ir + and K + yields from the quantum statistical value caused by 
the use of Boltzmann statistics as a function of the index k (h, I) (see Eq. [2]). 

Fig. [1] shows relative deviations of pion and kaon multiplicities from their quantum statis- 
tics values with increasing numbers of terms k in the expansion (J2]) . The calculations were 
performed for T = 157 MeV and V = 32 fm 3 , values that are relevant to freezeout condi- 
tions in e + e~ annihilation. It is clear from this figure that the Boltzmann approximation 
is by far not sufficient to reproduce the quantum statistics results. The pion yield under 
Boltzmann approximation deviates by more than 7% from the exact quantum statistics 
result. For kaons this difference is only 1%. While for kaons such an error is comparable 
to the 1 a error of the data in e + e~ collisions, for pions the deviation is significantly larger 
than the error in the data [26] . This underlines the importance of using quantum statistics 
for the calculation of multiplicities of light mesons. For pions, several terms are needed in 
the expansion ffTB"]) to achieve a precision well below 1%. It is also clear from Fig. [I] that 
the deviations depend on the mass of the particles and decrease quickly with increasing 
mass. For protons, the lightest fermions, the corresponding deviation from Fermi-Dirac 
statistics is below 0.1% already with only the first term, substantiating the applicability 
of the Boltzmann approximation for all fermions. 

In Fig. [2] we illustrate the relative error of the calculated multiplicity of pions and kaons 
using the partition function in the Boltzmann approximation as a function of VT 3 . De- 
viations from the exact quantum statistics values are seen to increase with VT 3 . This is 
due to the contribution of higher order terms in the expansion ffTBl) . which are suppressed 
for small VT 3 . For light bosons the largest deviations from the quantum statistics appear 
for large values of VT 3 i.e. when the system approaches GC thermodynamics. 
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Figure 2. Relative deviations of the pion and kaon multiplicities calculated using Boltzmann 
statistics from their quantum statistics values as a function of VT 3 . 

4 Conclusions 

We have presented a method to calculate the canonical partition function for the hadron 
resonance gas that accounts for exact conservation of baryon number, electric change, 
strangeness, charmness and bottomness. We have taken into account quantum statis- 
tics for light bosons and applied the Boltzmann approximation for fermions and heavier 
bosons. The results obtained here are an extension of previous studies which were re- 
stricted to the conservation of only three quantum numbers within the Boltzmann ap- 
proximation [24J. Our analytical expression for the partition function, which is represented 
as a series of Bessel functions, is stable in numerical implementations. It can be used for 
any value of the initial quantum numbers of the system and for arbitrary thermal pa- 
rameters. As an application of our results we have discussed the importance of quantum 
statistics in the calculations of pion and kaon multiplicities. The canonical partition func- 
tion which we have derived can be used to analyze different thermodynamical properties 
of the hadron resonance gas with the constraint of exact quantum number conservation. It 
can be also used to describe particle production in elementary and in heavy ion collisions 
within the statistical thermal model. 
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